Evaluating movement-based methods for estimating the frequency and timing of parturition in mule deer.

BACKGROUND: Information on reproduction of harvested species such as mule deer (Odocoileus hemionus) is vital for conservation and management. Furthermore, parturition in ungulates may be detected using patterns of movement logged by GPS transmitters. Several movement-based methods have been developed to detect parturition in ungulates including the Peterson method, behavioral change point analysis (BCPA), rolling minimum convex polygons (rMCP), individual-based method (IBM), and population-based method (PBM). Our objectives were to (1) test the accuracy and the precision of each previously described method and (2) develop an improved method optimized for mule deer that incorporated aspects of the other methods. METHODS: We determined parturition timing and status for female mule deer fitted with GPS collars and implanted with vaginal implant transmitters (VITs). We used movement patterns before and after parturition to set movement thresholds for each movement-based method. Following model training, we used location and birth date data from an external dataset to test the effectiveness of each movement-based method. Additionally, we developed a novel method for detecting parturition called the analysis of parturition indicators (API). We used two regression analyses to determine the accuracy and precision of estimates generated by each method. RESULTS: The six methods we employed varied in accuracy, with the API, rMCP, and BCPA being most accurate. Precision also varied among methods, with the API, rMCP, and PBM generating the most precise estimates of parturition dates. The API and the rMCP performed similarly and better overall than any of the other existing methods. CONCLUSIONS: We found that movement-based methods could be used to accurately and precisely detect parturition in mule deer. Further, we determined that the API and rMCP methods had the greatest overall success at detecting parturition in mule deer. The relative success of the API and rMCP may be attributed to the fact that both methods use home range size to detect parturition and are validated using known parturition dates of collared deer. We present the API as an efficient method of estimating birth status and timing of parturition of mule deer fitted with GPS transmitters, as well as affirm the effectiveness of a previously developed method, rMCP.

Estimating age of mule deer in the field: Can we move beyond broad age categories?

Age of individuals is an intrinsic demographic parameter used in the modeling and management of wildlife. Although analysis of cementum annuli from teeth is currently the most accurate method used to age ungulates, the age of live ungulates in the field can be estimated by examining tooth wear and tooth replacement patterns. However, there may be limitations to aging based on tooth wear as the rate of tooth wear likely varies among individuals due to factors such as age, diet, environment, and sex. Our objective was to determine the reliability of estimating age for mule deer based on tooth wear and tooth replacement patterns. We compared ages estimated by tooth wear (collected at time of capture for a statewide monitoring effort) to ages determined from cementum analysis (from teeth collected after mortalities of radio-tracked animals from the monitoring effort). Accuracy was high; ages estimated from tooth wear were within one year of cementum ages >75% of the time when aged by experienced observers. Bias in accuracy for estimates of age was low but slightly biased toward underestimation (i.e., 0.6 years on average)-especially as cementum age increased. Our results indicate that aging mule deer using patterns in tooth wear can be reliable if observers estimating age have experience using this method.

Variation in movement patterns of mule deer: have we oversimplified migration?

BACKGROUND: Conservation and management of migratory animals has gained attention in recent years, but the majority of research has focused on stereotypical 'migrant' and 'resident' behaviors, often failing to incorporate any atypical behaviors or characterize migratory behaviors beyond distance and timing of the migration. With migration threatened by anthropogenic development and climate change, it is crucial that we understand the full range of migratory behaviors. Our objective was to demonstrate and characterize the variation in migration strategies, including typical and atypical migratory behaviors for mule deer (Odocoileus hemionus) in Utah, USA. METHODS: Because calculation of common metrics such as distance, timing, and use of stopovers during migration did not adequately describe the variation we observed in migratory behavior for this species-particularly when animals visited multiple (> 3) ranges for extended lengths of time-we developed additional methods and categories to describe observed variation in migratory behavior. We first categorized trajectories based on the number of discrete, separate ranges and range shifts between them. Then, we further characterized the variation in migration strategies by examining the timing, duration, and distance traveled within each of the categories. We also examined if and how frequently individual deer switched among categories from year to year. RESULTS: We classified 1218 movement trajectories from 722 adult female mule deer, and found that 54.4% were dual-range migrants, who made one round-trip to one distinct range. Multi-range migrants (23.6%) made one round-trip during which they stayed at multiple discrete ranges. Commuters (1.0%) traveled to the same range multiple times, and poly migrants (1.5%) made multiple round-trips to different ranges. Gradual movers (2.5%) did not show a discrete range shift but moved gradually between ranges, whereas residents (12.6%) never left their home ranges, and dispersers (4.4%) left but never returned. Of the deer that we monitored for multiple years, 51.2% switched among categories. CONCLUSION: We conclude that the substantial number of atypical migratory strategies, as well as the number of deer that switched categories, underlines the importance of studying these less-stereotyped behaviors that may be exhibited by large proportions of populations. Acknowledging and investigating the full complexity and diversity in migratory strategies might uncover unknowns with respect to underlying factors and drivers of migration, and can help shape effective conservation strategies.

Spatiotemporal variation in the fecal microbiota of mule deer is associated with proximate and future measures of host health.

BACKGROUND: Mule deer rely on fat and protein stored prior to the winter season as an energy source during the winter months when other food sources are sparse. Since associated microorganisms ('microbiota') play a significant role in nutrient metabolism of their hosts, we predicted that variation in the microbiota might be associated with nutrient storage and overwintering in mule deer populations. To test this hypothesis we performed a 16S rRNA marker gene survey of fecal samples from two deer populations in the western United States before and after onset of winter. RESULTS: PERMANOVA analysis revealed the deer microbiota varied interactively with geography and season. Further, using metadata collected at the time of sampling, we were able to identify different fecal bacterial taxa that could potentially act as bioindicators of mule deer health outcomes. First, we identified the abundance of Collinsella (family: Coriobacteriaceae) reads as a possible predictor of poor overwintering outcomes for deer herds in multiple locations. Second, we showed that reads assigned to the Bacteroides and Mollicutes Order RF39 were both positively correlated with deer protein levels, leading to the idea that these sequences might be useful in predicting mule deer protein storage. CONCLUSIONS: These analyses confirm that variation in the microbiota is associated with season-dependent health outcomes in mule deer, which may have useful implications for herd management strategies.

De novo chromosome-length assembly of the mule deer (Odocoileus hemionus) genome.

The mule deer (Odocoileus hemionus) is an ungulate species that is distributed in a range from western Canada to central Mexico. Mule deer are an essential source of food for many predators, are relatively abundant, and commonly make broad migration movements. A clearer understanding of the mule deer genome can improve our knowledge of its population genetics, movements, and demographic history, aiding in conservation efforts. Their large population size, continuous distribution, and diversity of habitat make mule deer excellent candidates for population genomics studies; however, few genomic resources are currently available for this species. Here, we sequence and assemble the mule deer genome into a highly contiguous chromosome-length assembly for use in future research using long-read sequencing and Hi-C technologies. We also provide a genome annotation and compare demographic histories of the mule deer and white-tailed deer using the pairwise sequentially Markovian coalescent model. We expect this assembly to be a valuable resource in the continued study and conservation of mule deer.

The influence of habitat use on harvest vulnerability of cow elk (Cervus canadensis).

Pressure from hunting can alter the behavior and habitat selection of game species. During hunting periods, cervids such as elk (Cervus canadensis) typically select for areas further from roads and closer to tree cover, while altering the timing of their daily activities to avoid hunters. Our objective was to determine the habitat characteristics most influential in predicting harvest risk of elk. We captured 373 female elk between January 2015 and March 2017 in the Uinta-Wasatch-Cache National Forest and surrounding area of central Utah, USA. We determined habitat selection during the hunting season using a resource selection function (RSF) for 255 adult cow elk. Additionally, we used a generalized linear mixed model to evaluate risk of harvest based on habitat use within home ranges (3rd order selection) as well as the location of the home range on the landscape to evaluate vulnerability on a broader scale. Female elk selected for areas that reduced hunter access (rugged terrain, within tree cover, on private land). Age, elevation and distance to roads within a home range were most influential in predicting harvest risk (top model accounted for 36.2% of AIC weight). Elevation and distance to trees were most influential in predicting risk when evaluating the location of the home range (top model accounted for 42.1% of AIC weight). Vulnerability to harvest was associated with proximity to roads. Additionally, survival in our landscape decreased with age of femaleelk.

Fire, rodent herbivory, and plant competition: implications for invasion and altered fire regimes in the Mojave Desert.

Biological invasions are responsive to changing wildfire regimes related to human activities that are altering biological communities. Our objective was to investigate how fire, rodent activity, and competition among plant species modify plant community structure, invasion patterns, and vulnerability to altered fire regimes. We imposed experimental fires, and reduced rodent density using fencing in a full factorial design and quantified competitive interactions among plant species in the northeast Mojave Desert that has experienced dramatic increases in plant invasion and fire in recent years. Vegetation surveys were conducted in the experimental plots to determine plant density, cover, and biomass of herbaceous plants over a 5-year period. Rodent exclusion increased the density, cover, and biomass of Bromus rubens, an invasive annual grass, and density of forb species. In contrast, rodent exclusion decreased the density, cover, and biomass of Schismus spp. another dominant annual invader. Fire increased Schismus spp. and forb species density, cover, and biomass but decreased B. rubens density. Negative spatial correlation between B. rubens and Schismus spp., and forbs indicated interspecific competition among the dominant plant species. Fire reduced rodent community diversity (Shannon's) 2.5-fold, which was correlated with increases in B. rubens cover and biomass, and native forb diversity. Fire, high rodent diversity, and competition from the other plant species may decrease fire potential in our study system by reducing the density and biomass of B. rubens, which because of its taller growth form tends to ignite and carry fire better than Schismus spp. and forbs.

Spatial processes decouple management from objectives in a heterogeneous landscape: predator control as a case study.

Predator control is often implemented with the intent of disrupting top-down regulation in sensitive prey populations. However, ambiguity surrounding the efficacy of predator management, as well as the strength of top-down effects of predators in general, is often exacerbated by the spatially implicit analytical approaches used in assessing data with explicit spatial structure. Here, we highlight the importance of considering spatial context in the case of a predator control study in south-central Utah. We assessed the spatial match between aerial removal risk in coyotes (Canis latrans) and mule deer (Odocoileus hemionus) resource selection during parturition using a spatially explicit, multi-level Bayesian model. With our model, we were able to evaluate spatial congruence between management action (i.e., coyote removal) and objective (i.e., parturient deer site selection) at two distinct scales: the level of the management unit and the individual coyote removal. In the case of the former, our results indicated substantial spatial heterogeneity in expected congruence between removal risk and parturient deer site selection across large areas, and is a reflection of logistical constraints acting on the management strategy and differences in space use between the two species. At the level of the individual removal, we demonstrated that the potential management benefits of a removed coyote were highly variable across all individuals removed and in many cases, spatially distinct from parturient deer resource selection. Our methods and results provide a means of evaluating where we might anticipate an impact of predator control, while emphasizing the need to weight individual removals based on spatial proximity to management objectives in any assessment of large-scale predator control. Although we highlight the importance of spatial context in assessments of predator control strategy, we believe our methods are readily generalizable in any management or large-scale experimental framework where spatial context is likely an important driver of outcomes.

The influence of periodic increases of human activity on crepuscular and nocturnal mammals: Testing the weekend effect.

Human recreation can negatively affect wildlife, particularly on weekends when human activity is highest (i.e., the weekend effect). Much of what we understand about the weekend effect is based on research conducted on diurnal species, which have greater temporal overlap with humans. Because nocturnal species generally avoid times when humans are active, they are likely less affected by anthropogenic activity on weekends. Our objective was to test the weekend effect in relation to the degree of nocturnality of mammals in a recreational area. We predicted that as nocturnality increased, the effect of human activity would decrease. To address our objective, we placed 50 remote cameras along the Diamond Fork River in Utah from January to June 2015. We found that three out of the four focal species supported our predictions. Mule deer (crepuscular) reduced activity throughout our entire study area during weekends and avoided campgrounds. Beavers and mountain lions (both nocturnal) did not negatively respond to increased human activity. Raccoons (nocturnal) reduced activity during weekends, but only within campground areas. Our findings indicate that as the temporal overlap increases between wildlife and humans, so does the influence that humans have on wildlife.

A comparison of the effects of fire on rodent abundance and diversity in the Great Basin and Mojave Deserts.

As invasive grasses and fire increase in frequency and extent in North American deserts, they have the potential to affect animal communities through bottom-up forces. We experimentally tested the effects of fire on rodent communities of the Great Basin and Mojave Deserts. Fire decreased the abundance, richness, and diversity of rodents in the Great Basin after fire. In the Mojave, abundance was unaffected and diversity and species richness were greater on burned than unburned plots 4 months after fire. The effects of fire on rodent communities tended to decrease over time. The differences in effects between the deserts may be due to differences in the foraging preferences of the dominant species at each site. As these species are primarily herbivorous, short-term changes to the rodent community could have long-term implications by affecting the recovery of the plant community after fire.

Rodent herbivory differentially affects mortality rates of 14 native plant species with contrasting life history and growth form traits.

Ecosystems are transformed by changes in disturbance regimes including wildfire and herbivory. Rodent consumers can have strong top-down effects on plant community assembly through seed predation, but their impacts on post-germination seedling establishment via seedling herbivory need better characterization, particularly in deserts. To test the legacy effects of fire history, and native rodent consumers on seedling establishment, we evaluated factorial combinations of experimental exclusion of rodents and fire history (burned vs. unburned) on seedling survival of 14 native plant species that vary in their life history strategies and growth form in the Mojave Desert. Seedlings were placed into the experimental plots, and seedling survival was monitored daily for 8 days. The legacy effects of fire history had minimal effects on seedling survival, but rodent exclusion, year, and their interaction were strongly significant. Seedling survival rates were nearly sixfold greater in rodent exclusion plots compared to control plots in 2012 (53 vs. 9%) and 17-fold greater in 2013 (17 vs. 1%). The dramatic increase in seedling mortality from 2012 to 2013 was likely driven by an increase in rodent abundance and an outbreak of grasshoppers that appears to have intensified the rodent effect. There was strong variability in plant species survival in response to rodent herbivory with annual plants and forb species showing lower survival than perennial plants and shrub species. These results indicate that rodent consumers can strongly regulate seedling survival of native plant species with potentially strong regulatory effects on plant community development.

Increased Primary Production from an Exotic Invader Does Not Subsidize Native Rodents.

Invasive plants have tremendous potential to enrich native food webs by subsidizing net primary productivity. Here, we explored how a potential food subsidy, seeds produced by the aggressive invader cheatgrass (Bromus tectorum), is utilized by an important guild of native consumers--granivorous small mammals--in the Great Basin Desert, USA. In a series of field experiments we examined 1) how cheatgrass invasion affects the density and biomass of seed rain at the ecosystem-level; 2) how seed resources from cheatgrass numerically affect granivorous small mammals; and 3) how the food preferences of native granivores might mediate the trophic integration of cheatgrass seeds. Relative to native productivity, cheatgrass invasion increased the density and biomass of seed rain by over 2000% (P < 0.01) and 3500% (P < 0.01), respectively. However, granivorous small mammals in native communities showed no positive response in abundance, richness, or diversity to experimental additions of cheatgrass seeds over one year. This lack of response correlated with a distinct preference for seeds from native grasses over seeds from cheatgrass. Our experiments demonstrate that increased primary productivity associated with exotic plant invasions may not necessarily subsidize consumers at higher trophic levels. In this context, cheatgrass invasion could disrupt native food webs by providing less-preferred resources that fail to enrich higher trophic levels.

Negative effects of an exotic grass invasion on small-mammal communities.

Exotic invasive species can directly and indirectly influence natural ecological communities. Cheatgrass (Bromus tectorum) is non-native to the western United States and has invaded large areas of the Great Basin. Changes to the structure and composition of plant communities invaded by cheatgrass likely have effects at higher trophic levels. As a keystone guild in North American deserts, granivorous small mammals drive and maintain plant diversity. Our objective was to assess potential effects of invasion by cheatgrass on small-mammal communities. We sampled small-mammal and plant communities at 70 sites (Great Basin, Utah). We assessed abundance and diversity of the small-mammal community, diversity of the plant community, and the percentage of cheatgrass cover and shrub species. Abundance and diversity of the small-mammal community decreased with increasing abundance of cheatgrass. Similarly, cover of cheatgrass remained a significant predictor of small-mammal abundance even after accounting for the loss of the shrub layer and plant diversity, suggesting that there are direct and indirect effects of cheatgrass. The change in the small-mammal communities associated with invasion of cheatgrass likely has effects through higher and lower trophic levels and has the potential to cause major changes in ecosystem structure and function.

Water developments and canids in two North American deserts: a test of the indirect effect of water hypothesis.

Anthropogenic modifications to landscapes intended to benefit wildlife may negatively influence wildlife communities. Anthropogenic provisioning of free water (water developments) to enhance abundance and distribution of wildlife is a common management practice in arid regions where water is limiting. Despite the long-term and widespread use of water developments, little is known about how they influence native species. Water developments may negatively influence arid-adapted species (e.g., kit fox, Vulpes macrotis) by enabling water-dependent competitors (e.g., coyote, Canis latrans) to expand distribution in arid landscapes (i.e., indirect effect of water hypothesis). We tested the two predictions of the indirect effect of water hypothesis (i.e., coyotes will visit areas with free water more frequently and kit foxes will spatially and temporally avoid coyotes) and evaluated relative use of free water by canids in the Great Basin and Mojave Deserts from 2010 to 2012. We established scent stations in areas with (wet) and without (dry) free water and monitored visitation by canids to these sites and visitation to water sources using infrared-triggered cameras. There was no difference in the proportions of visits to scent stations in wet or dry areas by coyotes or kit foxes at either study area. We did not detect spatial (no negative correlation between visits to scent stations) or temporal (no difference between times when stations were visited) segregation between coyotes and kit foxes. Visitation to water sources was not different for coyotes between study areas, but kit foxes visited water sources more in Mojave than Great Basin. Our results did not support the indirect effect of water hypothesis in the Great Basin or Mojave Deserts for these two canids.

Long-lasting effects of maternal condition in free-ranging cervids.

Causes of phenotypic variation are fundamental to evolutionary ecology because they influence the traits acted upon by natural selection. One such cause of phenotypic variation is a maternal effect, which is the influence of the environment experienced by a female (and her corresponding phenotype) on the phenotype of her offspring (independent of the offspring's genotype). While maternal effects are well documented, the longevity and fitness impact of these effects remains unclear because it is difficult to follow free-living individuals through their reproductive lifetimes. For long-lived species, it has been suggested that maternal effects are masked by environmental variables acting on offspring in years following the period of dependence. Our objective was to use indirect measures of maternal condition to determine if maternal effects have long-lasting influences on male offspring in two species of cervid. Because antlers are sexually selected, we used measures of antler size at time of death, 1.5-21.5 years after gestation to investigate maternal effects. We quantified antler size of 11,000 male elk and mule deer born throughout the intermountain western US (6 states) over nearly 30 years. Maternal condition during development was estimated indirectly using a suite of abiotic variables known to influence condition of cervids (i.e., winter severity, spring and summer temperature, and spring and summer precipitation). Antler size of male cervids was significantly associated with our indirect measure of maternal condition during gestation and lactation. Assuming the correctness of our indirect measure, our findings demonstrate that antler size is a sexually selected trait that is influenced-into adulthood-by maternal condition. This link emphasizes the importance of considering inherited environmental effects when interpreting population dynamics or examining reproductive success of long-lived organisms.

Vigilance of kit foxes at water sources: a test of competing hypotheses for a solitary carnivore subject to predation.

Animals that are potential prey do not respond equally to direct and indirect cues related to risk of predation. Based on differential responses to cues, three hypotheses have been proposed to explain spatial variation in vigilance behavior. The predator-vigilance hypothesis proposes that prey increase vigilance where there is evidence of predators. The visibility-vigilance hypothesis suggests that prey increase vigilance where visibility is obstructed. Alternatively, the refuge-vigilance hypothesis proposes that prey may perceive areas with low visibility (greater cover) as refuges and decrease vigilance. We evaluated support for these hypotheses using the kit fox (Vulpes macrotis), a solitary carnivore subject to intraguild predation, as a model. From 2010 to 2012, we used infrared-triggered cameras to record video of kit fox behavior at water sources in Utah, USA. The refuge-vigilance hypothesis explained more variation in vigilance behavior of kit foxes than the other two hypotheses (AICc model weight=0.37). Kit foxes were less vigilant at water sources with low overhead cover (refuge) obstructing visibility. Based on our results, the predator-vigilance and visibility-vigilance hypotheses may not be applicable to all species of prey. Solitary prey, unlike gregarious prey, may use areas with concealing cover to maximize resource acquisition and minimize vigilance.